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Abstract: Recently we have studied the Lorentzian version of the IIB matrix model 
as a nonperturbative formulation of superstring theory. By Monte Carlo simulation, we 
have shown that the notion of time — as well as space — emerges dynamically from this 
model, and that we can uniquely extract the real-time dynamics, which turned out to be 
rather surprising: after some "critical time", the S0(9) rotational symmetry of the nine- 
dimensional space is spontaneously broken down to SO (3) and the three-dimensional space 
starts to expand rapidly. In this paper, we study the same model based on the classical 
equations of motion, which are expected to be valid at later times. After providing a general 
prescription to solve the equations, we examine a class of solutions, which correspond to 
manifestly commutative space. In particular, we find a solution with an expanding behavior 
that naturally solves the cosmological constant problem. 
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1. Introduction 

There are many fundamental questions in cosmology, which can, in principle, be answered 
by superstring theory. Describing the birth of the universe is one of the most fundamental 
ones. It is recognized, however, that the cosmic singularity is not resolved generally in 
perturbative string theory [1-5]. Therefore, in order to study the very early universe, 
we definitely need a nonperturbative formulation. Among various proposals [6-8] based 
on matrix models,^ the IIB matrix model [7] looks most natural for describing the birth 
of the universe, since not only space but also time is expected to appear dynamically 
from the matrix degrees of freedom. Also the model is unique in that it is a manifestly 
covariant formulation, while the other proposals are based on the light-cone formulation, 
which breaks covariance. 

^For earlier attempts to apply matrix models to cosmology, see ref. [9-19]. 
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One of the important issues in the IIB matrix model is to identify the configurations 
of matrices that dominate the path integral and to determine the corresponding space- 
time structure. This was studied by various approaches [20-34] in the Euclidean version of 
the model, which was shown to have finite partition function without any cutoffs [35,36]. 
However, the Euclidean model does not seem suitable for cosmology since it does not 
provide the real-time dynamics. Furthermore, a recent study based on the gaussian expan- 
sion method suggests that the space-time obtained dynamically in the Euclidean model is 
three-dimensional rather than four-dimensional [34] . While this conclusion itself may have 
a profound implication, we do not know yet how it should be physically interpreted. 

All these considerations led us to study the Lorentzian version of the IIB matrix model 
nonperturbatively [37]. By Monte Carlo simulation, wc have shown that the Lorentzian 
model can be made well-defined nonperturbatively by first introducing infrared cutoffs 
and then removing them appropriately in the large- limit. We have also found that the 
eigenvalue distribution of the matrix in the temporal direction extends in that limit, which 
implies that time emerges dynamically. (Supersymmetry of the model plays a crucial role.) 
Indeed we were able to extract a unique real-time dynamics, which turns out to have a 
surprising property. After some critical time, the S0(9) rotational symmetry of the space 
is broken spontaneously down to S0(3), and the three-dimensional space starts to expand 
rapidly. This result can be interpreted as the birth of the universe. Note that the concept 
of time-evolution also emerges dynamically in this model without ever having to specify 
the initial condition. 

Although the length of time-evolution that has been extracted from Monte Carlo sim- 
ulation is restricted due to finite N, we consider that the whole history of the universe 
can be obtained from the same model in the large- A limit. If this is true, we should be 
able to answer various important questions in both cosmology and particle physics. For 
instance, we may obtain the microscopic description of inflation. If we can reach the time 
at which stringy excitations and quantum gravitational effects become negligible, we may 
see how particles in the Standard Model (or possibly its extension) starts to appear. If we 
can study the behaviors of the model at much later times, wc may be able to understand 
why the expansion of our universe is accelerating in the present epoch. Finally we may 
even predict how our universe will be like in the future. 

While the late-time behaviors are difficult to study by direct Monte Carlo methods, 
the classical equations of motion are expected to become more and more valid at later 
times since the value of the action increases with the cosmic expansion. We will see that 
there are actually many classical solutions, which is reminiscent of the fact that string 
theory possesses infinitely many vacua that are perturbatively stable. However, unlike 
in perturbative string theory, we have the possibility to pick up the unique solution that 
describes our universe by requiring smooth connection to the behavior at earlier times 
accessible by Monte Carlo simulation. From this perspective, we consider it important to 
classify the classical solutions and to examine their cosmological implications. The aim of 
this paper is to make a first step in that direction. In particular, we find a classical solution 
with an expanding behavior that can naturally solve the cosmological constant problem. 

Another issue we would like to address in this paper concerns how the commutative 
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space-time appears from this model. This is important since the S0(9) symmetry breaking 
observed in the Monte Carlo simulation is understood intuitively by a mechanism, which 
relies crucially on the fact that the matrices that represent the space-time are noncommu- 
tative [37]. We show that classical solutions which correspond to manifestly commutative 
space can be easily constructed, and discuss the vanishing of noncommutativity between 
space and time in some simple examples. This implies that the emergence of commutative 
space-time is indeed possible in this model at later times. 

The rest of this paper is organized as follows. In section 2 we briefly review the IIB 
matrix model and the results obtained by our previous Monte Carlo studies [37] . In section 
3 we provide a general prescription to find classical solutions. In particular, we show 
that this can be done systematically by using Lie algebras. Then we restrict ourselves 
to manifestly space-space commutative solutions and present a complete classification of 
such solutions within certain simplifying ansatz. The particular solution discussed in our 
previous publication [38] also appears in this classification. In section 4 we obtain explicit 
time-evolution of the scale factor and the Hubble parameter for some simple solutions, and 
discuss their cosmological implications. Section 5 is devoted to a summary and discussions. 
In appendix A we review the irreducible unitary representations of the SU(1, 1) algebra, 
which will be needed in constructing the solutions discussed in section 4. In appendix B 
we list the other simple solutions which are not discussed in section 4. In appendix C we 
give some examples of solutions that are not manifestly space-space commutative. 

2. The IIB matrix model and the birth of the universe 

In this section we discuss why the IIB matrix model is considered as a nonperturbative 
formulation of type IIB superstring theory in ten dimensions, and how the results suggesting 
the birth of the universe were obtained by our previous Monte Carlo studies of this model. 
The action of the IIB matrix model is given as [7] 



where An {/j, = 0, • • • ,9) and (cc = 1, ■ " " > 16) are N x N Hermitian matrices. The 

Lorentz indices fi and v are contracted using the metric r] = diag(— 1, 1, • • • ,1). The 16 x 16 
matrices are ten-dimensional gamma matrices after the Weyl projection, and the unitary 
matrix C is the charge conjugation matrix. The action has manifest SO (9,1) symmetry, 
where and ^q, transform as a vector and a Majorana-Weyl spinor, respectively. 

There are various evidences that the model gives a nonperturbative formulation of 
superstring theory. First of all, the action (2.1) can be viewed as a matrix regularization 
of the worldshcet action of type IIB superstring theory in a particular gauge known as the 
Schild gauge [7]. It has also been argued that configurations of block-diagonal matrices 
correspond to a collection of disconnected worldsheets with arbitrary genus. Therefore, 
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instead of being equivalent just to the worldsheet theory, the large-A^ limit of the matrix 
model is expected to be a second-quantized theory of type IIB superstrings, which includes 
multi-string states. Secondly, D-branes are represented as classical solutions in the matrix 
model, and the interaction between them calculated at one loop reproduced correctly the 
known results from type IIB superstring theory [7] . Thirdly, one can derive the light-cone 
string field theory for the type IIB case from the matrix model [39] with a few assumptions. 
In the matrix model, one can define the Wilson loops, which can be naturally identified 
with the creation and annihilation operators of strings. Then, from the Schwinger-Dyson 
equations for the Wilson loops, one can actually obtain the string field Hamiltonian. 

In all these connections to string theory, it is crucial that the model has two kinds of 
fermionic symmetries given as 

f 5^^) A, = ie^CT,^ , 

i ''"'^^ = « ' (2 5) 

=621, 

where 1 is the unit matrix. It also has the bosonic symmetry given by 

[sm =0. ^ 

Let us denote the generators of (2.4), (2.5) and (2.6) by Q^^\ Q^'^^ and P^, respectively, 
and define their linear combinations 

Q(i)=Qa)+Q(2) , Q(2)=,(Q{i)_g(2)). (2.7) 
Then, we find that the generators satisfy the algebra 

[eiCQ(^), e2CQ(^)] = -2(5^^eiCr'^e2P^ , (2.8) 

where i, j = 1, 2. This is nothing but the ten-dimensional J\f = 2 supersymmetry. It is 
known that field theories with this symmetry necessarily include gravity, which suggests 
that so docs the IIB matrix model. When we identify (2.8) with the ten-dimensional 
M = 2 supersymmetry, the symmetry (2.6) is identified with the translational symmetry 
in ten dimensions, which implies that the eigenvalues of A^j^ should be identified with 
the coordinates of ten-dimensional space-time. This identification is consistent with the 
one adopted in stating the evidences listed in the previous paragraph, and shall be used 
throughout this paper as well. 

An interesting feature of the IIB matrix model is that the space-time itself is treated 
as a part of dynamical degrees of freedom in the matrices. One can therefore try to 
identify the dominant configurations of matrices in the path integral and to determine the 
corresponding space-time structure. For that purpose, one needs to define the partition 
function as a finite matrix integral. This is nontrivial since the bosonic part of the action 
(2.2) is not bounded from below. By decomposing it into two terms 

= - 2^ tr {Foif + ^ tr {F,^f , (2.9) 
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where we have dGfincd Hcriniticin nicitrices = i[A^,A^], we find that the first term is 
negative, whereas the second term is positive. A common way to overcome this problem 
is to make a Wick rotation, which amounts to making replacements Aq i-> iAio, r° i-> 
— zr^°. The Euclidean version of the model obtained in this way is well-defined nonpertur- 
batively without any cutoffs [35,36]. The space-time structure has been studied by various 
approaches in this Euclidean model [20-34]. However, it is nontrivial whether the Wick 
rotation is valid in a theory including gravity [40,41]. Furthermore, the recent result based 
on the gaussian expansion method [34] suggests that the space-time appearing dynamically 
in the Euclidean model is three-dimensional rather than four-dimensional. 

For these reasons, we studied the Lorentzian model nonperturbatively for the first time 
in ref . [37] . In order to make the partition function finite, we introduced the infrared cutoffs 
by imposing the following constraints 



where k and L are the cutoff parameters. We have shown by Monte Carlo simulation 
that these cutoffs can be removed in the large- iV limit in such a way that the physical 
quantities scale. The resulting theory thus obtained has no parameters other than one 
scale parameter. This feature is precisely what one expects for nonperturbative string 
theory. 

It turned out that not only space but also time emerges dynamically in this Lorentzian 
model. We found that the eigenvalue distribution of Aq extends in the large-AT limit. 
Here supersymmetry of the model plays a crucial role. If we omit fermions, the eigenvalue 
distribution has a finite extent. Furthermore, it turned out that one can extract the real- 
time dynamics by working in the SU(A^) basis which diagonalizes Aq. We found that after 
a critical time, the SO (9) symmetry is spontaneously broken down to the SO (3), and three 
out of nine spatial directions start to expand. This can be interpreted as the birth of the 
universe. Note that the real-time dynamics is an emergent notion in this model, and we 
do not even have to specify the initial conditions. The above result is unique in that sense. 

3. Classical solutions in the Lorentzian model 
3.1 General prescription to find classical solutions 

In this subsection we present a general prescription to find classical solutions in the 
Lorentzian model. It is important to take into account that the two cutoffs had to be 
introduced in order to make the model well-defined as we reviewed in the previous section. 
Since the inequalities (2.10) are actually saturated as is also seen by Monte Carlo simula- 
tion [37], we search for stationary points of the bosonic action S}^ for fixed ;^tr (^o)^ and 
;^tr(^i)2. Then we have to extremize the function 



^tr(^o)' < i^j^tiiAif , 
ltr(A)2 < l2 , 



(2.10) 




(3.1) 
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where A and A are the Lagrange multiphcrs. 

Differentiating (3.1) with respect to Aq and Ai, we obtain^ 

-[Ao, [Ao, Ai]] + [Aj, [Aj,Ai]] - A^ = , (3.2) 
[Aj,[Aj,Ao]]-~XAo = , (3.3) 

respectively. Differentiating (3.1) with respect to A and A, we obtain 

^tT{Al) = kL' , (3.4) 

^tr(Af ) = , (3.5) 

respectively. Once wc obtain a solution to eqs. (3.2) and (3.3), we can substitute them into 
(3.4) and (3.5) to determine A and A as a function of k and L. 

Here we point out that the terms in eqs. (3.2) and (3.3) proportional to the Lagrange 
multipliers break the S0(9, 1) symmetry in general. In fact, we will see that there is a 
set of solutions with A = A, which do not suffer from this explicit breaking. However, we 
do not impose this condition from the outset in order to keep our analysis as general as 
possible. 

In what follows, we consider solutions with 

Ai = for z > d , (3.6) 

where d < 9. This is motivated from our observation in Monte Carlo simulation that Ai 
in the extra dimensions remain small when the three-dimensional space expands. Prom 
this point of view, one may think that we should choose d = 3. However, one can actually 
construct a solution with larger d by taking a direct sum of solutions with smaller d, say 
with d = 1, as we will see in 3.3.2. We therefore keep d arbitrary for the moment. 

A general prescription to solve the equations of motion (3.2) and (3.3) within the 
ansatz (3.6) is given as follows. Let us first define a sequence of commutation relations 



[Ai,Aj]=iCij , (3.7) 

[A,, C,k] = iDijk , (3.8) 

[Ao,Ai]=iEi, (3.9) 

[Ao,E,]=iFi, (3.10) 

[A,Ej]=iGij , ••• , (3.11) 



where 1 < i,j,k < d and the symbols on the right-hand side represent Hermitian operators 
newly defined. Then we determine the relationship among -^i, Cij, Dijk, Ei, Fi, Gij, ■ ■ ■ 
so that the equations of motion (3.2) and (3.3) and the Jacobi identities are satisfied. We 
obtain a Lie algebra in this way. Considering that all the operators are Hermitian, each 
unitary representation of the Lie algebra gives a classical solution. 

^Classical solutions in the IIB matrix model have been studied in ref. [42] for A = A = 0. They are 
studied in the Euclidean version in ref. [43]. 
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3.2 Lie algebra for manifestly commutative space 

Since we expect that commutative space appears at later times, we restrict ourselves to 
solutions corresponding to manifestly commutative space in what follows. (See appendix 
C for examples of solutions without this restriction.) This implies that we impose 

[Ai,Aj] = 0. (3.12) 

Then we follow the general prescription described in section 3.1. In particular, we show 
that one can actually close the algebra with a finite number of generators by imposing a 
simple condition. 

Let us consider the relationship among Aq, Ai, Cij, Dijk, Ei, Fi, Gij in (3.7) ~ (3.11). 
First of aU, (3.12) implies that 

Ci,=0, A,fc = 0. (3.13) 

It is convenient to make the irreducible decomposition of a d-dimensional second-rank 
tensor Gij in (3.11) as 

= M,j + + ^6ijH , (3.14) 



where = Mji , ^ Mjj = , (3.15) 

i=l 

Nij = -Nji . (3.16) 

Prom the equations of motion (3.2) and (3.3), we obtain 

Fi = XAi , (3.17) 
H = \Ao , (3.18) 

where we have used (3.9), (3.10), (3.11), (3.12) and (3.14). 
Next we consider the Jacobi identities 

[Ao, [Ai,Aj]] + [A,, [Aj,Ao]] + [Aj, [Ao,Ai]] = , (3.19) 

[^0, [Ai,Ej]] + [A, [Ej,Ao]] + [Ej, [Ao,A,]] = , (3.20) 

[Ai, [Ej,Ek]] + [E,, [Ek, A,]] + [Ek, [A„ E^]] = , (3.21) 

[E,, [A,,Ak]] + [Aj, [Ak,E,]] + [Ak, [Ei,A,]] = . (3.22) 

Using (3.9), (3.11), (3.12) and (3.14) in (3.19), we find 

Nij = . (3.23) 

Using (3.17) and (3.18) in (3.20), we find 

[Ao,Mij] = 0, [Ei,Ej] = 0. (3.24) 
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Similarly, from (3.21) and (3.22), we obtain 

[Ej^Mki] - i^SkiAj = [Ek, Mij] - i^SijAk , (3.25) 

[Aj,Mik] - i^6kiEj = [Ak, Mij] - i^SijEk , (3.26) 
respectively. We can easily verify that the Jacobi identities 

[Ao, [Ei, Ej]] + [Ei, [Ej,Ao]] + [Ej, [Ao,Ei]] = , (3.27) 

[Ai, [Aj,Ak]] + [Aj, [Ak,Ai]] + [Ak, [A, Aj]] = , (3.28) 

[Ei, [Ej,Ek]] + [Ej, [Ek, Ei]] + [Ek, [Ei, Ej]] = (3.29) 

are trivially satisfied, hence giving no new relations among the operators. 

Now that all the Jacobi identities among Aq, Ai and Ei are satisfied, let us move on to 
the Jacobi identities that include Mij. In general we need to introduce some new operators, 
which appear from the commutator of Mij and one of Aq, Ai, Ei. One way to close the 
algebra without introducing new operators is to impose that Mij is diagonal. Let us denote 
the diagonal elements as Mj = Ma, which satisfy 

d 

J2Mi = (3.30) 

1=1 

due to the traceless condition (3.15). In what follows, we will see that Aq, Ai, Ei and Mj 
form a Lie algebra. 

First, we find that (3.25) and (3.26) lead to 

[Ei, Mj] = - ddij)Ai , (3.31) 

[Ai,Mj]=i^{l-d6ij)Ei . (3.32) 
Applying these commutation relations to the Jacobi identity 

[Mi, [Aj,Ek]] + [Aj, [Ek, Mi]] + [Ek, [Mi, Aj]] = , (3.33) 

we obtain 

[Mi,Mj]=0. (3.34) 
To summarize, the commutation relations among Aq, Ai, Ei and Mj are obtained as 

[Ai, Aj] = , [Ao, Ai] = iEi , [Aq, Ei] = iXAi , 

[Ei,Ej] = 0, [Ai,Ej]=i6ij(^^Ao + Mi^ , [Ao,Mi]=0, 

[Ai, Mj] = ^^(1 - d6ij)Ei , [Ei, Mj] = i^(l - dSij)Ai , [Mi, Mj] = , (3.35) 

where 1 < i, j < d with 1 < d < 9. It is straightforward to verify that all the remaining 
Jacobi identities including Mj are satisfied due to the commutation relations (3.35). Thus 
we obtain the Lie algebra (3.35), which gives a class of manifestly space-space commutative 
solutions. 
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3.3 Some simplifications of the Lie algebra 

In this subsection we consider some special cases of (3.35), which corresponds to simple 
Lie algebras. 

3.3.1 the case with Mj = and A = 

First we point out that one can set Mj = and A = consistently in eq. (3.35) for arbitrary 
d, which results in the Lie algebra 

[Ai, Aj] = , [Ao, Ai] = iE, , [^0, Ei] = iXAi, 

[Ei,Ej\=Q, [Ai,Ej\=0. (3.36) 

We can further simplify (3.36) by setting E^ = ±-\/A^i, which yields^ 

[Ai, Aj] = , [Ao, Ai] = ±^^^XAi . (3.37) 

The classical solutions obtained from this Lie algebra were studied in ref . [38] . 

3.3.2 the d = \ case 

The Lie algebra simplifies considerably in the d = 1 case. In this case, eq. (3.30) implies 
that Ml = 0. Thus eq. (3.35) reduces to 

[Ao, Ai] = iE , [Aq, E] = iXAi , [Ai,E] = fXAo , (3.38) 

where we define E = Ei. Note that the A = case of (3.38) is identical to the d = 1 case 
of (3.36). 

As we mentioned in section 3.1, we can use the d = 1 solutions (3.38) to construct 
new solutions representing a higher-dimensional space-time in the following way. For that, 
we note that the equations of motion (3.2) and (3.3) have S0(9) symmetry. Rotating the 
solution (3.38) by an S0(9) transformation, we obtain an equivalent solution, which has 
the i-th spatial matrix given by riAi {i = 1, ■ ■ ■ ,9) with rf = 1. Taking a direct sum of 
these solutions with various values of , we obtain a new solution: 

A'o = ^0 <^ Ifc , 

4 = Ai C5 diag(rf ^ , rf ) , • • • , rf ^) , (3.39) 
where rj™^^ = 1 {m = 1, ■ ■ ■ ,k) , (3.40) 

and Ife is the k x k unit matrix. In particular, we can construct an SO{D) symmetric 
solution by requiring that r^™) 's be distributed uniformly on a unit S"^-!, where l< D <9. 
The D = 4 case would then be a physically interesting solution which represent (3 -|- 1)- 
dimensional space-time with R x geometry. 

The Lie algebra (3.37) in the d = 3 and d = 4 cases correspond to Ai% and A^'^f in Table I of ref. [44], 

respectively. 
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Let us consider the case in which A 7^ and A 7^ 0. In this case, the Lie algebra (3.38) 
can be identified either with the SU(1, 1) algebra (the SL(2, R) algebra) 

[To, Ti] = iT2 , [T2, To] = in , [Ti, T2] = -in , (3.41) 
or with the SU(2) algebra 

[Li, L2] = iL3 , [L2, Ls] = zLi , [L3, Li] = iL2 , (3.42) 
depending on the signs of A and A as follows.^ 

(a) A > and A > : SU(1, 1) algebra 

^0 = aT2 , Ai = bTo, E = cTi , 

A = , A = 6^ , ab = c. (3.43) 

(b) A < and A < : SU(1, 1) algebra 

Ao = aTo , Ai = bTi, E = cTs , 

A = -a^ , A = -6^ , ab = c. (3.44) 

(c) A > and A < : SU(1, 1) algebra 

^0 = aT2 , Ai = bTi , E = cTo , 

A = , A = -6^ , ab = c. (3.45) 

(d) A < and A > : SU(2) algebra 

^0 = , Ai = bLi , E = CL2 , 

A = -a^ , A = 6^ , ab = c. (3.46) 

The cases in which A = or A = are discussed in appendix B. Thus we find that the 
solutions with d = 1 are classified into the nine cases; namely, (a)~(d) in this subsection 
and (i)~(v) in appendix B. In section 4 we discuss cosmological implications of the above 
four cases (a)~(d). 

3.3.3 the d = 2 case 

As a more complicated example, we discuss the d = 2 case of (3.35). Since this case will 
not be discussed further in this paper, impatient readers may jump into section 4. 

Let us note that eq. (3.30) gives Mi = —M2 = M. For A 7^ and A 7^ 0, we can rescale 
Aq, Ai, Ei and M appropriately so that eq. (3.35) can be rewritten as 

[^1,^2] = [E^,E2] = [Ai,E2] = [A2,Ei] = [Ao,M] = , 



*These two algebras appear in Table I of ref. [44] as A^^s and ^3,9, respectively. 
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[^0, ^i] = iEi , [vlo, ^2] = iE2 , [Ao, Ei] = isign(A)^i , [Aq, E2] = isign(A)^2 , 
[Ai, El] = 2i(sign(A)^o + M) , [^2,^2] = 2i(sign(A)^o - M) , 
[Ai,M] = -isign(A)Si , [A2, M] = isignC\)E2 , 

[Ei,M] = -isign(AA)^i , [E2, M] = isign(AA)A2 , (3.47) 

where sign( • ) represents the sign function. We compare the above algebra with the 
S0(2, 2) algebra and S0(4) algebra, which can be expressed in a unified way as 

[Lai3, L^s] = iga'yLjSS + ^5/35-^07 " iQasLfi^ - igp^L^s , (3.48) 
where a, /3, 7, (5 = 1,2,3,4 and ga/3 represents the metric. By making an identification 

^0 = L23 , 

Ai = L12 + -L34 , A2 = L12 — -L34 , 

El = 522 -£'13 — 533-f'24 , E2 = 522-^13 + 533-^^24 , 

M = -sign(A)Li4 , (3.49) 



it is easy to check that eq. (3.47) is satisfied if 



(gap) 



diag(l, 1,-1,-1) forA>0, A>0 

diag(l, -1, -1, 1) for A < 0, A < 

diag(l, -1, 1,-1) for A > 0, A < 

^ diag(l, 1, 1, 1) for A < 0, A > 



(3.50) 



The first three cases correspond to the SO (2, 2) algebra, whereas the last one corresponds 
to the S0(4) algebra. A unitary representation of these algebras corresponds to a classical 
solution. More in-depth studies of these solutions are left for future investigations. 



4. Cosmological implications of some simple classical solutions 

In this section we discuss the cosmological implications of the SU(1, 1) solutions (a)~(c) 
and the SU(2) solution (d) discussed in section 3.3.2. As a warming up, we start with the 
SU(2) solution, which is simpler, and then move on to the SU(1, 1) solutions, which exhibit 
physically more interesting behaviors. 

4.1 Solutions based on the SU(2) algebra — a warm-up 

Let us consider the solution (3.46) based on the SU(2) algebra (3.42). As we mentioned 
in section 3.1, we use the unitary representations. The irreducible unitary representations 
of the SU(2) algebra is specified by their spins J, which are non-negative integers or half- 
integers. In the spin J representation, the matrix elements of the generators in eq. (3.42) 
are given by 

iLi)rnn = ^ V( J - n) ( J + n + l)5m,n+l + ^\/(J + n)(J - n + l)Sm,n-l , 
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{L2)mn = ^-V {J - 'n'){J + n + l)Sm,n+l - (J + n){J -n+ l)5m,n-l , 



(4.1) 



where — J < m,n < J. From eq. (3.46), we find that Aq is diagonal, whereas Ai has a 
tri-diagonal structure. This motivates us to extract the time evolution of space from the 
3x3 submatrices of and Ai defined as [37] 



Ao{n) = a 



/n- 1 \ 

n 
y n + ly 

/ ^y{J + n){J-n + l) \ 

y(J + n)(J-n + l) y(J-n)(J + n+l) 

^/{J -n){J + n + \) 



(4.2) 



(4.3) 



where — J + l<n<J — 1. 

Similarly, we consider an S0(4) symmetric solution (3.39), where rf^^'^ = 1 (m 
1, • • • ,k) and r^"*) are uniformly distributed on a unit S^. In this case, we define 



A'^in) = Ao{n) ® Ife , 
A'i{n) = ^(n) (g)diag(rf\' 



(4.4) 



where ^^(n) represents the structure of space at a discrete time n. We also define the 
extent of space at a discrete time n by 



1 



R(n)^\l-triA[in))^ 



62 



(J(J + l)-n2) . 



(4.5) 



Let us then discuss the continuum limit, in which wc send a ^ and J oo, and 
define the continuum time hj t = na. We see from eq. (4.5) that a nontrivial dependence 
of on f is obtained by keeping tmax = Jet and | = a fixed. This leads to 



R{t) = i?maxW 1 



where = 



atr. 



(4.6) 



From the fact that Aq is diagonal and Ai has a tri-diagonal structure, we consider 
that the space-time noncommutativity vanishes in the continuum limit. Let us also note 
that the extents of space and time defined by (3.4) and (3.5) are given by L ~ aimax and 
K ~ 1/a. 

Let us next discuss the cosmological implication of this solution. For that we naively 
identify R{t) with the scale factor in the Friedman-Robertson- Walker metric. Then the 
space-time is i? x S"^, where the radius of S"^ is given by R{t). Figure 1 shows the time 
dependence of R{t). The universe expands towards t = and shrinks after t = 0. The 
Hubble parameter can be defined in terms of the scale factor R{t) as 



Hit) 



m 
m 



cR{t) 



-1(1+"') 



(4.7) 
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where c is a constant. From this, we can evaluate the parameter w using 

Let us recall that w = 1/3, w = and w = —1 correspond to the radiation dominated 
universe, the matter dominated universe and the cosmological constant term, respectively. 
In the present case of SU(2) solution, we find from (4.6) that 

H = ^^^Rl^-R^, (4.9) 



i-max 



w = %^-i (4.10) 

for t < 0. The parameter w is w = 1/3 dX t = —tmax, and it diverges as one approaches 
i = 0. 




Figure 1: The time dependence of R{t) in the SU(2) solution 



4.2 Solutions based on the SU(1, 1) algebra 

Next we discuss the solutions (a), (b) and (c) in section 3.3.2, which are based on the 
SU(1,1) algebra (3.41). We construct S0(4) symmetric solutions for these cases as in 
eq. (3.39). The irreducible unitary representations of the SU(1, 1) algebra are summarized 
in appendix A. Apart from the trivial representation, there are three types: the primary 
unitary series representation (PUSR) (A. 4), the discrete series representation (DSR) (A. 6) 
and (A. 7), the complementary unitary series representation (CUSR) (A. 5). In all these 
representations, the generator Tq is diagonal. Hence, we can analytically treat the case 
(b), in which is proportional to Tq, while we need to diagonalize Aq numerically in the 
cases (a) and (c). For this reason, we start with the case (b). In this case, we find from 
(A. 3) that Ai has a tri-diagonal structure. Therefore, we extract 3x3 submatrices AqIu) 
and Ai{n) similarly to the SU(2) case discussed in section 4.1. In what follows, we discuss 
the solutions for each representation separately. 
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Let us first discuss the solutions corresponding to PUSR. (n) and Ai takes the form 



fn-l+e 



Ao{n) = a 








\ 

n + e 
n + l + ej 



Aiin) 



ib 
2" 







— n + ip + 2 — e 




n + ip — ^ + e \ 

n + ip + ^ + e 





—n + ip — ly — e 



(4.11) 



(4.12) 



where e = or ^ and p is a non-negative number, which specifies a representation. 

When we consider an S0(4) symmetric sohition, we define Ag(n) and A'j^{n) as in 
eq. (4.4). Then we find that the extent of space R{n) at a discrete time n becomes 



(4.13) 



Let us take the continuum limit. We define the continuum time hy t = na and take 
the limit in which a — > with | = a fixed. We can tune p so that to = pa is fixed. Then 
R{t) is given by 



R{t) 



(4.14) 



As in the SU(2) case, the continuum limit of this solution represents a commutative 
(3+l)-dimensional space-time. In order to evaluate L and k in (3.4) and (3.5), we introduce 
a cutoff N for the dimension of the representation. Then, it is easy to see that (3.4) and 
(3.5) give L ~ Na and k ~ 1/a. The infrared cutoff L is removed by sending N to infinity 
faster than 1/a. The other parameter k, which corresponds to the ratio of tr(^o)^/-^ 
to ti {Ai)'^ /N , is finite in the continuum limit, which looks different from the situation 
encountered in Monte Carlo studies, where we had to send k to infinity [37]. This is not 
so surprising, however, given that we are looking at different time regions, and the speed 
of expansion changes qualitatively depending on which region we are looking at. 

Here we naively identify R{t) with the scale factor again. Then the space-time is 
R X S^, where the radius R{t) of is time dependent. Prom eq. (4.14), we obtain the 
Hubble parameter H and the parameter w as 



H 



a 



w = — 



2f2 



i?2 



3t2 



1 
3 



(4.15) 
(4.16) 



We find that w converges to — | as t — >■ oo, which corresponds to the expansion of universe 
with a constant velocity. 

If we identify Iq with the present time, the present value of w is —1. This value of w 
corresponds to the cosmological constant, which explains the present accelerating expansion 
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of the universe. Moreover, the corresponding cosmological constant becomes of the order of 
(1/io)^, which naturally solves the cosmological constant problem. As we mentioned above, 
w increases with time and approaches — |. This means that the cosmological constant 
actually vanishes in the future. In fig. 2 we show the time dependence of R{t) and the 
parameter w. 




Figure 2: The time dependence of the scale factor R{t) (Left) and the parameter w (Right) in the 
SU(1, 1) solution with the PUSR. 



Let us next discuss the solutions corresponding to DSR. Ao{n) takes the form (4.11), 
where e = or ^ , and Ai takes the form 



ib 



v^(n + r)(n-r-l) 

Mn) = ^ I -^y{n + T){n-T -1) ^(n - r)(n + r + 1) 

-y(n-r)(n + T + l) 



(4.17) 



with r = —1, —2, —3, • • • for e = and r = — ^, — |, — |, • • • for e = ^. There is a constraint 
n < — r — e + 1 or n<r — e — 1. Then R{n) is given by 



R{n) = ^jin^-Tir + 1)). (4.18) 
In the continuum limit, we can tune r so that to = arij + 1) is fixed. Then R{t) is given 

by 



R{t) = ^-{t^-tl), (4.19) 

where the range of t is restricted to either i > io or i < —Iq. 

Prom eq. (4.19), we obtain the Hubble parameter H and the parameter w as 



H=^^\ R^ + ^ , (4.20) 
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w 



0^2 
3t2 



(4.21) 



for t > tf). We find that the parameter w becomes ttJ = 1/3 at t = to ^ind = at t = \/2io- 
The former corresponds to the radiation dominant universe, while the latter to the matter 
dominant universe. Thus this solution may represent some part of the history of the 
universe. Figure 3 shows the time dependence of R{t) and the parameter w, respectively. 



1/3 




Figure 3: The time dependence of the scale factor R{t) (Left) and the parameter w (Right) in the 
SU(1, 1) solution with the DSR. 



Next we discuss the solutions corresponding to CUSR. AqIh) takes the form (4.11), 
where e = 0, and Ai{n) takes the form (4.17), where r within — 1 < r < specifies a 
representation. The extent of space R{n) is given by eq. (4.18). In the continuum limit, 
we obtain (4.19), where to = 0- 

Finally, we discuss the cases (a) and (c) in section 3.3.2. In these cases, we have 
diagonalized Aq numerically. The physical consequence obtained from (a) is essentially the 
same as (b). On the other hand, no accelerating expansion is obtained for the case (c). 
R{t) has two peaks, one in the t < region and the other in the t > region, and the 
minimum lies at t = 0. 

5. Summciry and discussions 

In this paper we studied the late time behaviors of the universe in the Lorentzian version 
of the IIB matrix model. We investigated the classical equations of motion, which are 
expected to be valid at later times. This is a complementary approach to Monte Carlo 
simulation,^ which was used previously to study the birth of the universe in the same 
model. First we provided a general prescription to solve the equations of motion. The 
problem reduces to that of finding a unitary representation of a Lie algebra. 

^See ref. [45] for a recent review on Monte Carlo studies of matrix models and supersymmetric gauge 
theories in the context of string theory. 
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In this way, we obtained a class of solutions that are manifestly space-space commu- 
tative. The simplest ones in this class are the d = 1 solutions with A2 = ■ ■ ■ = Aq = 0, 
from which we can easily construct the ones representing higher dimensional space-time 
as well. We made a complete classification of such solutions. Some solutions represent 
expanding (3 -|- l)-dimensional universe without space-time noncommutativity in the con- 
tinuum limit. In particular, we find that there exists a solution, in which the parameter 
w changes smoothly from —1 to —1/3. This explains why we seem to have a tiny cos- 
mological constant in the present epoch, and hence can naturally solve the cosmological 
constant problem. While we do not insist that this particular solution really describes our 
universe, we consider that the cosmological constant problem can be naturally solved in 
the Lorentzian matrix model in a similar manner. 

Corresponding to what we have done in Monte Carlo simulation, wc have introduced 
infrared cutoffs in both the temporal and spatial directions. These are represented by the 
Lagrange multipliers A and A introduced in the action (3.1). In general, this breaks the 
S0(9,l) symmetry of the model explicitly. Let us note, however, that it is possible to 
have A = A in the cases (a) and (b) in section 3.3.2. Such solutions break the SO (9,1) 
symmetry spontaneously. We expect that the explicit breaking of the Lorentz symmetry 
by the infrared cutoffs disappears in the large- limit. If that is really the case, we should 
select a solution with A = A. It is intriguing to note that the cases (a) and (b) are indeed 
the ones that are physically interesting. 

The (3-|-l)-dimensional space-time represented by the solutions discussed mainly in 
this paper has the topology R x S^. This is a restriction which we have as long as we 
construct such solutions based on the d = 1 solution. In other constructions, we can also 
obtain solutions representing a space with the topology of a three-dimensional ball as we 
discussed in appendix B. While the Monte Carlo results seem to be more consistent with 
the latter topology of the space, it remains to be seen what kind of topology is actually 
realized at later times. 

Below we list some directions for future investigations. 

First we consider it important to examine the stability of the solutions we found in 
this paper. It would be also interesting to calculate the one-loop effective action around 
the solutions. That would tell us the validity of the solutions, and we should be able to 
know how late the time should be for the solutions to be valid. 

Secondly it is important to understand better how one should extract the information 
of the space-time metric from a matrix configuration. Ref. [46] shows that this is indeed 
possible, in principle, if one interprets the matrix as a covariant derivative on the space- 
time manifold, where the general coordinate invariance is realized manifestly as a subgroup 
of the SU(A^) symmetry. However, this interpretation is different from the one adopted 
in this paper, which is compatible with the supersymmetry as we reviewed in section 2. 
The precise relationship between the two interpretations is yet to be clarified, although 
it is tempting to consider that they are related to each other by T-duality of type IIB 
supcrstring theory. In this work wc have naively identified R{t) with the scale factor in the 
Friedman-Robertson- Walker metric when we discuss cosmological implications in section 
4. It remains to be seen whether this identification can somehow be justified. 
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Thirdly we consider it important to study a wider class of solutions using the general 
prescription provided in this paper. In particular, it would be interesting to examine the 
solutions, which are not manifestly space-space commutative, based on the Lie algebra 
(3.49) or the one given in appendix C. Also it would be interesting to investigate solutions 
with nontrivial structure in the extra dimensions. Such structure is expected to play a 
crucial role [47, 48] in determining the matter content at late times and in finding how 
the standard model appears from the matrix model. Eventually, wc have to single out the 
solution, which is smoothly connected to the unique result at earlier times accessible by 
Monte Carlo simulation. 

Developments in the above directions would enable us to solve various fundamental 
problems in particle physics and cosmology. For instance, we should be able to understand 
the mechanism of inflation and to clarify what the dark matter and the dark energy are. 
We hope that the present work will trigger such developments. 
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A. Uniteiry representations of the SU(1, 1) algebra 

In this appendix we summarize the unitary representations of the SU(1, 1) algebra (3.41) 
based on ref. [49]. The generators are realized in the space of square integrable functions 
in the region [0, 27r], which we denote as £^(0, 27r) in what follows. They are given as 



(T + e)e^^ + (r-e)e-^^-2sin^-^ 

da ^ 

-(r + e)e^^ + (r - e)e-^^ - 2i cos 6^ 

do 



(A.1) 



where < ^ < 27r, and r G C and e G M are parameters. It is easy to verify that these 
operators satisfy the SU(1, 1) algebra 

[%,Ti] = iT2, [T2,To]=iTi, [Ti,T2] = -iTo . (A.2) 

Taking the set of functions {e"*'"^; m G Z} as a basis of i2^(0, 27r), one obtains the matrix 
elements of the generators as 

iTo)mn = (e + n)5mn , 



-18- 



ijl)mn = --2in-T + e)5m,n+l + -jijl + T + e)5m,n-l , 

(7i)m„ = -^{n - T + e)6m,n+l " ^{n + T + e)5m,n-l , (A.3) 

where m,n E 1^. 

The unitary irreducible representations are classified as follows. We denote the matrix 
elements of T^u in (3.41) by (T]u)mn5 which differs from (7^)mrn in general, due to some 
factors introduced to define the scalar product that realizes the unitarity. Since SU(1, 1) 
is a noncompact group, all the nontrivial unitary representations are infinite dimensional. 

1) primary unitary series representations 

T = ip- - (p G M>o) , e = or - , 

2) complementary unitary series representations 

-1<T<1, e=0. 



r(T — Tfi ~\~ i)r( — T — ?T/) 

{Tp)mn = I TTTnT T J 0~iJ,)mn ■ (A. 5) 

\r{T — n + l)T{—T — m) ' 



3) discrete series representations (I) 



r = -l,-2,-3,--- , e = , 

or 

13 5 1 



"2 ' 2' 2'"' ' ^ 2 ' 



_ f T{T + m + e + l)T{-T + n + e) Y 
^ l^r(r + n + e + l)r(-r + m + e)y' ' 

m,n> —T — e . (A.6) 

4) discrete series representations (II) 

r = -l,-2,-3,--- , e = 0, 
or 

__1 3 5 _1 
^~ ~2'~2'~2''" ' ^~ 2 ' 



_ ^ r(T-m-e + l)r(-r-n-e) V,^, 
^ Vr(T-n-e+l)r(-r-m-e)y' ' 

m,n<r-e. (A.7) 

5) trivial representation 

= . (A.8) 
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B. The other classical solutions in the d = 1 case 

In this appendix we discuss the solutions based on the Lie algebra (3.38) obtained for d = 1 
when A = or A = 0. The solutions can be classified into the following five cases. ^ 

An interesting feature of the A = case is that wc can construct solutions representing 
higher dimensional space-time with topology other than S^^^. The reason is that we can 
rescale Ai and E in (3.38) without changing A when A = 0. Therefore, we do not need to 
impose the condition (3.40) in constructing the new solution (3.39). For instance, we can 
distribute r^-"^^ uniformly in a three-dimensional ball so that the solution represents a 
(3 + l)-dimensional space-time with SO (3) symmetry. 

(i) A = and A = 

The nontrivial irreducible representations are parametrized by a G M (a 7^ 0). ^0 and 
Ai are given by the operators acting on the space of functions of x with integrability, 
which is denoted by >C^(R) in what follows. The operators are given explicitly as 

Ao = —iaVX-^ , Ai = X , E = —a . (B.l) 

ax 

This solution represents an infinitely long static D-string. 

(ii) A > and A = 

In the nontrivial irreducible representations, Aq, Ai and E are operators acting on >C^(M), 
which are given by 

Ao = —iVX-^ , yli = a cosh a:; -I- 6 sinh X , = — '\/A(asinha; -|- 6cosha;) , (B.2) 
ax 

where a, 6 G M. 

For a = b, the solution reduces to Ai = aexp{x) and E = VXAi, which corresponds 
to the d = 1 case of (3.37). Wc can construct an S0(3) symmetric solution from this case 
by the aforementioned procedure. The solution thus obtained is equivalent to the one we 
constructed from the Lie algebra (3.37) with d = 3 in ref. [38], which represents a (3 + 1)- 
dimensional expanding universe. For a / 5, we can also construct an SO (3) symmetric 
solution in the same way. We have checked numerically that the resulting solution exhibits 
essentially the same behavior as the one for a = b. 

(ni) A < and A = 

In the nontrivial irreducible representations, ^1 and E are operators acting on C^{0, 27r), 
which are given as 

Ao = —iV—X-^ , Ai = acosx + bsinx , E = V— A(a sin x — 6 cos x) . (B.3) 



®The five cases (i)~(v) correspond to ^3,1, ^3,4, ^43,6, ^4.3, 6 and ^3,4, respectively, in Table I of ref. [44]. 
In particular, ^43,1 is a nilpotent Lie algebra. 
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This case is analytically tractable. For the SO (3) symmetric solution obtained from the 
Lie algebra (B.3), we can define R{t) as in the case of the SU(2) and SU(1, 1) solutions. 
We find that R{t) = constant. 

(iv) A = and A > 

In the nontrivial irreducible representations, Ao, Ai and E are operators acting on £^{0, 27r), 
which are given as 

Aq = acosx + bsinx , Ai = —iV^^ , E = \/\(—asmx + bcosx) . (B.4) 

ax 

We have seen numerically that the S0(3) symmetric solution obtained from (B.4) exhibits 
almost the same behavior as the one in (iii). 



(v) A = and A < 

In the nontrivial irreducible representations, Aq, Ai and E are operators acting on L^(R) : 



= a cosh a:; + 6 sinh X , Ai = —iy —X— , E = y —X{asmh.x + bcoshx) . (B.5) 



We have seen numerically that R{t) in the S0(3) symmetric solution obtained from (B.5) 
exhibits a behavior different from (ii). The solutions for 6 = has an expanding regime 
only, whereas the solutions for a = have both expanding and contracting regimes. 

C. Examples of classical solutions describing noncommutative space 

In this appendix we present some examples of classical solutions which are not manifestly 
space-space commutative. These solutions are based on the Lie algebras S0(6), S0(5, 1) 
and S0(4, 2), and we interpret them as describing (3 + l)-dimensional universes with S0(4) 
symmetry. These Lie algebras obey the commutation relations 





(C.l) 



where a, /S, 7, 5 = 1, 2, • • • , 



6. The non-vanishing components of ga/s are 




gii = l (i = 1,2,3,4), 

_ / 1 for S0(6), S0(5,l) , 



(C.2) 



We set 




4) , 



(C.3) 
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Then it is easy to verify that the equations (3.2) and (3.3) are satisfied if 



A = -a 555566 + 46 5-55 , 

A = 46^555 . (C.4) 
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